Wolbachia infection at least partially rescues the fertility and ovary defects of several new Drosophila melanogaster bag of marbles protein-coding mutants.

The D. melanogaster protein coding gene bag of marbles (bam) plays a key role in early male and female reproduction by forming complexes with partner proteins to promote differentiation in gametogenesis. Like another germline gene, Sex lethal, bam genetically interacts with the endosymbiont Wolbachia, as Wolbachia rescues the reduced fertility of a bam hypomorphic mutant. Here, we explored the specificity of the bam-Wolbachia interaction by generating 22 new bam mutants, with ten mutants displaying fertility defects. Nine of these mutants trend towards rescue by the wMel Wolbachia variant, with eight statistically significant at the fertility and/or cytological level. In some cases, fertility was increased a striking 20-fold. There is no specificity between the rescue and the known binding regions of bam, suggesting wMel does not interact with one singular bam partner to rescue the reproductive phenotype. We further tested if wMel interacts with bam in a non-specific way, by increasing bam transcript levels or acting upstream in germline stem cells. A fertility assessment of a bam RNAi knockdown mutant reveals that wMel rescue is specific to functionally mutant bam alleles and we find no obvious evidence of wMel interaction with germline stem cells in bam mutants.

Wolbachia infection at least partially rescues the fertility and ovary defects of several new Drosophila melanogaster bag of marbles protein-coding mutants.

The D. melanogaster protein coding gene bag of marbles ( bam ) plays a key role in early male and female reproduction by forming complexes with partner proteins to promote differentiation in gametogenesis. Like another germline gene, Sex lethal , bam genetically interacts with the endosymbiont Wolbachia , as Wolbachia rescues the reduced fertility of a bam hypomorphic mutant. Here, we explored the specificity of the bam-Wolbachia interaction by generating 22 new bam mutants, with ten mutants displaying fertility defects. Nine of these mutants trend towards rescue by the w Mel Wolbachia variant, with eight statistically significant at the fertility and/or cytological level. In some cases, fertility was increased a striking 20-fold. There is no specificity between the rescue and the known binding regions of bam , suggesting w Mel does not interact with one singular bam partner to rescue the reproductive phenotype. We further tested if w Mel interacts with bam in a non-specific way, by increasing bam transcript levels or acting upstream in germline stem cells. A fertility assessment of a bam RNAi knockdown mutant reveals that w Mel rescue is specific to functionally mutant bam alleles and we find no obvious evidence of w Mel interaction with germline stem cells in bam mutants. Author Summary: Reproduction in the Drosophila melanogaster fruit fly is dependent on the bag of marbles ( bam ) gene, which acts early in the process of generating eggs and sperm. Mutations to this gene negatively impact the fertility of the fly, causing it to be sterile or have fewer progeny. Interestingly, we find that the bacteria Wolbachia , which resides within reproductive cells across a wide range of insects, partially restores the fertility and ovary phenotype of several bam mutants of which the resultant Bam protein is altered from wildtype. The protein function of Bam is further suggested to be important by the lack of rescue for a fly that has a fertility defect due to low expression of a non-mutated bam gene. Previous work makes similar conclusions about Wolbachia with another reproductive gene, Sex lethal ( Sxl ), highlighting the potential for rescue of fertility mutants to occur in a similar way across different genes. An understanding of the ways in which Wolbachia can affect host reproduction provides us with context with which to frame Wolbachia 's impact on host genes, such as bam and Sxl, and consider the evolutionary implications of Wolbachia 's infection in D. melanogaster fruit flies.

Evolution under a model of functionally buffered deleterious mutations can lead to positive selection in protein-coding genes.

Selective pressures on DNA sequences often result in departures from neutral evolution that can be captured by the McDonald-Kreitman (MK) test. However, the nature of such selective forces often remains unknown to experimentalists. Amino acid fixations driven by natural selection in protein-coding genes are commonly associated with a genetic arms race or changing biological purposes, leading to proteins with new functionality. Here, we evaluate the expectations of population genetic patterns under a buffering mechanism driving selective amino acids to fixation, which is motivated by an observed phenotypic rescue of otherwise deleterious nonsynonymous substitutions at bag of marbles (bam) and Sex lethal (Sxl) in Drosophila melanogaster. These two genes were shown to experience strong episodic bursts of natural selection potentially due to infections of the endosymbiotic bacteria Wolbachia observed among multiple Drosophila species. Using simulations to implement and evaluate the evolutionary dynamics of a Wolbachia buffering model, we demonstrate that selectively fixed amino acid replacements will occur, but that the proportion of adaptive amino acid fixations and the statistical power of the MK test to detect the departure from an equilibrium neutral model are both significantly lower than seen for an arms race/change-in-function model that favors proteins with diversified amino acids. We find that the observed selection pattern at bam in a natural population of D. melanogaster is more consistent with an arms race model than with the buffering model.

Wolbachia genetically interacts with the bag of marbles germline stem cell gene in male D. melanogaster.

The bacterial endosymbiont Wolbachia manipulates reproduction of its arthropod hosts to promote its own maternal vertical transmission. In female D. melanogaster , Wolbachia has been shown to genetically interact with three key reproductive genes ( bag of marbles ( bam ) , Sex-lethal, and mei-P26) , as it rescues the reduced female fertility or fecundity phenotype seen in partial loss-of-function mutants of these genes . Here, we show that Wolbachia also partially rescues male fertility in D. melanogaster carrying a new, largely sterile bam allele when in a bam null genetic background. This finding shows that the molecular mechanism of Wolbachia 's influence on its hosts' reproduction involves interaction with genes in males as well as females, at least in D. melanogaster .

Functional Divergence of the bag-of-marbles Gene in the Drosophila melanogaster Species Group.

In Drosophila melanogaster, a key germline stem cell (GSC) differentiation factor, bag of marbles (bam) shows rapid bursts of amino acid fixations between sibling species D. melanogaster and Drosophila simulans, but not in the outgroup species Drosophila ananassae. Here, we test the null hypothesis that bam's differentiation function is conserved between D. melanogaster and four additional Drosophila species in the melanogaster species group spanning approximately 30 million years of divergence. Surprisingly, we demonstrate that bam is not necessary for oogenesis or spermatogenesis in Drosophila teissieri nor is bam necessary for spermatogenesis in D. ananassae. Remarkably bam function may change on a relatively short time scale. We further report tests of neutral sequence evolution at bam in additional species of Drosophila and find a positive, but not perfect, correlation between evidence for positive selection at bam and its essential role in GSC regulation and fertility for both males and females. Further characterization of bam function in more divergent lineages will be necessary to distinguish between bam's critical gametogenesis role being newly derived in D. melanogaster, D. simulans, Drosophila yakuba, and D. ananassae females or it being basal to the genus and subsequently lost in numerous lineages.

Recommendations for improving statistical inference in population genomics.

The field of population genomics has grown rapidly in response to the recent advent of affordable, large-scale sequencing technologies. As opposed to the situation during the majority of the 20th century, in which the development of theoretical and statistical population genetic insights outpaced the generation of data to which they could be applied, genomic data are now being produced at a far greater rate than they can be meaningfully analyzed and interpreted. With this wealth of data has come a tendency to focus on fitting specific (and often rather idiosyncratic) models to data, at the expense of a careful exploration of the range of possible underlying evolutionary processes. For example, the approach of directly investigating models of adaptive evolution in each newly sequenced population or species often neglects the fact that a thorough characterization of ubiquitous nonadaptive processes is a prerequisite for accurate inference. We here describe the perils of these tendencies, present our consensus views on current best practices in population genomic data analysis, and highlight areas of statistical inference and theory that are in need of further attention. Thereby, we argue for the importance of defining a biologically relevant baseline model tuned to the details of each new analysis, of skepticism and scrutiny in interpreting model fitting results, and of carefully defining addressable hypotheses and underlying uncertainties.

Molecular population genetics of Sex-lethal (Sxl) in the Drosophila melanogaster species group: a locus that genetically interacts with Wolbachia pipientis in Drosophila melanogaster.

Sex-lethal (Sxl) is the sex determination switch in Drosophila, and also plays a critical role in germ-line stem cell daughter differentiation in Drosophila melanogaster. Three female-sterile alleles at Sxl in D. melanogaster were previously shown to genetically interact to varying degrees with the maternally inherited endosymbiont Wolbachia pipientis. Given this genetic interaction and W. pipientis' ability to manipulate reproduction in Drosophila, we carried out a careful study of both the population genetics (within four Drosophila species) and molecular evolutionary analysis (across 20 Drosophila species) of Sxl. Consistent with earlier studies, we find that selective constraint has played a prominent role in Sxl's molecular evolution within Drosophila, but we also observe patterns that suggest both episodic bursts of protein evolution and recent positive selection at Sxl. The episodic nature of Sxl's protein evolution is discussed in light of its genetic interaction with W. pipientis.

Diverse wMel variants of Wolbachia pipientis differentially rescue fertility and cytological defects of the bag of marbles partial loss of function mutation in Drosophila melanogaster.

In Drosophila melanogaster, the maternally inherited endosymbiont Wolbachia pipientis interacts with germline stem cell genes during oogenesis. One such gene, bag of marbles (bam) is the key switch for differentiation and also shows signals of adaptive evolution for protein diversification. These observations have led us to hypothesize that W. pipientis could be driving the adaptive evolution of bam for control of oogenesis. To test this hypothesis, we must understand the specificity of the genetic interaction between bam and W. pipientis. Previously, we documented that the W. pipientis variant, wMel, rescued the fertility of the bamBW hypomorphic mutant as a transheterozygote over a bam null. However, bamBW was generated more than 20 years ago in an uncontrolled genetic background and maintained over a balancer chromosome. Consequently, the chromosome carrying bamBW accumulated mutations that have prevented controlled experiments to further assess the interaction. Here, we used CRISPR/Cas9 to engineer the same single amino acid bam hypomorphic mutation (bamL255F) and a new bam null disruption mutation into the w1118 isogenic background. We assess the fertility of wildtype bam, bamL255F/bamnull hypomorphic, and bamL255F/bamL255F mutant females, each infected individually with 10 W. pipientis wMel variants representing three phylogenetic clades. Overall, we find that all of the W. pipientis variants tested here rescue bam hypomorphic fertility defects with wMelCS-like variants exhibiting the strongest rescue effects. In addition, these variants did not increase wildtype bam female fertility. Therefore, both bam and W. pipientis interact in genotype-specific ways to modulate female fertility, a critical fitness phenotype.

Baker's Yeast Clinical Isolates Provide a Model for How Pathogenic Yeasts Adapt to Stress.

Global outbreaks of drug-resistant fungi such as Candida auris are thought to be due at least in part to excessive use of antifungal drugs. Baker's yeast Saccharomyces cerevisiae has gained importance as an emerging opportunistic fungal pathogen that can cause infections in immunocompromised patients. Analyses of over 1000 S. cerevisiae isolates are providing rich resources to better understand how fungi can grow in human environments. A large percentage of clinical S. cerevisiae isolates are heterozygous across many nucleotide sites, and a significant proportion are of mixed ancestry and/or are aneuploid or polyploid. Such features potentially facilitate adaptation to new environments. These observations provide strong impetus for expanding genomic and molecular studies on clinical and wild isolates to understand the prevalence of genetic diversity and instability-generating mechanisms, and how they are selected for and maintained. Such work can also lead to the identification of new targets for antifungal drugs.

The importance of the Neutral Theory in 1968 and 50 years on: A response to Kern and Hahn 2018.

A recent article reassessing the Neutral Theory of Molecular Evolution claims that it is no longer as important as is widely believed. The authors argue that "the neutral theory was supported by unreliable theoretical and empirical evidence from the beginning, and that in light of modern, genome-scale data, we can firmly reject its universality." Claiming that "the neutral theory has been overwhelmingly rejected," they propose instead that natural selection is the major force shaping both between-species divergence and within-species variation. Although this is probably a minority view, it is important to evaluate such claims carefully in the context of current knowledge, as inaccuracies can sometimes morph into an accepted narrative for those not familiar with the underlying science. We here critically examine and ultimately reject Kern and Hahn's arguments and assessment, and instead propose that it is now abundantly clear that the foundational ideas presented five decades ago by Kimura and Ohta are indeed correct.

Incompatibilities in Mismatch Repair Genes MLH1-PMS1 Contribute to a Wide Range of Mutation Rates in Human Isolates of Baker's Yeast.

Laboratory baker's yeast strains bearing an incompatible combination of MLH1 and PMS1 mismatch repair alleles are mutators that can adapt more rapidly to stress, but do so at the cost of long-term fitness. We identified 18 baker's yeast isolates from 1011 surveyed that contain the incompatible MLH1-PMS1 genotype in a heterozygous state. Surprisingly, the incompatible combination from two human clinical heterozygous diploid isolates, YJS5845 and YJS5885, contain the exact MLH1 (S288c-derived) and PMS1 (SK1-derived) open reading frames originally shown to confer incompatibility. While these isolates were nonmutators, their meiotic spore clone progeny displayed mutation rates in a DNA slippage assay that varied over a 340-fold range. This range was 30-fold higher than observed between compatible and incompatible combinations of laboratory strains. Genotyping analysis indicated that MLH1-PMS1 incompatibility was the major driver of mutation rate in the isolates. The variation in the mutation rate of incompatible spore clones could be due to background suppressors and enhancers, as well as aneuploidy seen in the spore clones. Our data are consistent with the observed variance in mutation rate contributing to adaptation to stress conditions (e.g., in a human host) through the acquisition of beneficial mutations, with high mutation rates leading to long-term fitness costs that are buffered by mating or eliminated through natural selection.

Mismatch Repair Incompatibilities in Diverse Yeast Populations.

An elevated mutation rate can provide cells with a source of mutations to adapt to changing environments. We identified a negative epistatic interaction involving naturally occurring variants in the MLH1 and PMS1 mismatch repair (MMR) genes of Saccharomyces cerevisiae We hypothesized that this MMR incompatibility, created through mating between divergent S. cerevisiae, yields mutator progeny that can rapidly but transiently adapt to an environmental stress. Here we analyzed the MLH1 and PMS1 genes across 1010 S. cerevisiae natural isolates spanning a wide range of ecological sources (tree exudates, Drosophila, fruits, and various fermentation and clinical isolates) and geographical sources (Europe, America, Africa, and Asia). We identified one homozygous clinical isolate and 18 heterozygous isolates containing the incompatible MMR genotype. The MLH1-PMS1 gene combination isolated from the homozygous clinical isolate conferred a mutator phenotype when expressed in the S288c laboratory background. Using a novel reporter to measure mutation rates, we showed that the overall mutation rate in the homozygous incompatible background was similar to that seen in compatible strains, indicating the presence of suppressor mutations in the clinical isolate that lowered its mutation rate. This observation and the identification of 18 heterozygous isolates, which can lead to MMR incompatible genotypes in the offspring, are consistent with an elevated mutation rate rapidly but transiently facilitating adaptation. To avoid long-term fitness costs, the incompatibility is apparently buffered by mating or by acquiring suppressors. These observations highlight effective strategies in eukaryotes to avoid long-term fitness costs associated with elevated mutation rates.

Recent and Long-Term Selection Across Synonymous Sites in Drosophila ananassae.

In Drosophila, many studies have examined the short- or long-term evolution occurring across synonymous sites. Few, however, have examined both the recent and long-term evolution to gain a complete view of this selection. Here we have analyzed Drosophila ananassae DNA polymorphism and divergence data using several different methods, and have identified evidence of positive selection favoring preferred codons in both recent and long-term evolutionary time scale. Further in D. ananassae, the strength of selection for preferred codons was stronger on the X chromosome compared to the autosomes. We show that this stronger selection is not due to higher gene expression of X-linked genes. Analysis of the selectively neutral introns indicated that the X chromosome also had a preference for GC over AT nucleotides, potentially from GC-biased gene conversions (gcBGCs) that can also affect the base composition of synonymous sites. Thus selection for preferred codons and gcBGC both seem to be partially responsible for shaping the D. ananassae synonymous site evolution.

Inversions and adaptation to the plant toxin ouabain shape DNA sequence variation within and between chromosomal inversions of Drosophila subobscura.

Adaptation is defined as an evolutionary process allowing organisms to succeed in certain habitats or conditions. Chromosomal inversions have the potential to be key in the adaptation processes, since they can contribute to the maintenance of favoured combinations of adaptive alleles through reduced recombination between individuals carrying different inversions. We have analysed six genes (Pif1A, Abi, Sqd, Yrt, Atpα and Fmr1), located inside and outside three inversions of the O chromosome in European populations of Drosophila subobscura. Genetic differentiation was significant between inversions despite extensive recombination inside inverted regions, irrespective of gene distance to the inversion breakpoints. Surprisingly, the highest level of genetic differentiation between arrangements was found for the Atpα gene, which is located outside the O1 and O7 inversions. Two derived unrelated arrangements (O3+4+1 and O3+4+7) are nearly fixed for several amino acid substitutions at the Atpα gene that have been described to confer resistance in other species to the cardenolide ouabain, a plant toxin capable of blocking ATPases. Similarities in the Atpα variants, conferring ouabain resistance in both arrangements, may be the result of convergent substitution and be favoured in response to selective pressures presumably related to the presence of plants containing ouabain in the geographic locations where both inversions are present.

Molecular Evolution of Drosophila Germline Stem Cell and Neural Stem Cell Regulating Genes.

Here, we study the molecular evolution of a near complete set of genes that had functional evidence in the regulation of the Drosophila germline and neural stem cell. Some of these genes have previously been shown to be rapidly evolving by positive selection raising the possibility that stem cell genes as a group have elevated signatures of positive selection. Using recent Drosophila comparative genome sequences and population genomic sequences of Drosophila melanogaster, we have investigated both long- and short-term evolution occurring across these two different stem cell systems, and compared them with a carefully chosen random set of genes to represent the background rate of evolution. Our results showed an excess of genes with evidence of a recent selective sweep in both germline and neural stem cells in D. melanogaster. However compared with their control genes, both stem cell systems had no significant excess of genes with long-term recurrent positive selection in D. melanogaster, or across orthologous sequences from the melanogaster group. The evidence of long-term positive selection was limited to a subset of genes with specific functions in both the germline and neural stem cell system.

A Genetic Incompatibility Accelerates Adaptation in Yeast.

During mismatch repair (MMR) MSH proteins bind to mismatches that form as the result of DNA replication errors and recruit MLH factors such as Mlh1-Pms1 to initiate excision and repair steps. Previously, we identified a negative epistatic interaction involving naturally occurring polymorphisms in the MLH1 and PMS1 genes of baker's yeast. Here we hypothesize that a mutagenic state resulting from this negative epistatic interaction increases the likelihood of obtaining beneficial mutations that can promote adaptation to stress conditions. We tested this by stressing yeast strains bearing mutagenic (incompatible) and non-mutagenic (compatible) mismatch repair genotypes. Our data show that incompatible populations adapted more rapidly and without an apparent fitness cost to high salt stress. The fitness advantage of incompatible populations was rapid but disappeared over time. The fitness gains in both compatible and incompatible strains were due primarily to mutations in PMR1 that appeared earlier in incompatible evolving populations. These data demonstrate a rapid and reversible role (by mating) for genetic incompatibilities in accelerating adaptation in eukaryotes. They also provide an approach to link experimental studies to observational population genomics.

Population Genomics of Infectious and Integrated Wolbachia pipientis Genomes in Drosophila ananassae.

Coevolution between Drosophila and its endosymbiont Wolbachia pipientis has many intriguing aspects. For example, Drosophila ananassae hosts two forms of W. pipientis genomes: One being the infectious bacterial genome and the other integrated into the host nuclear genome. Here, we characterize the infectious and integrated genomes of W. pipientis infecting D. ananassae (wAna), by genome sequencing 15 strains of D. ananassae that have either the infectious or integrated wAna genomes. Results indicate evolutionarily stable maternal transmission for the infectious wAna genome suggesting a relatively long-term coevolution with its host. In contrast, the integrated wAna genome showed pseudogene-like characteristics accumulating many variants that are predicted to have deleterious effects if present in an infectious bacterial genome. Phylogenomic analysis of sequence variation together with genotyping by polymerase chain reaction of large structural variations indicated several wAna variants among the eight infectious wAna genomes. In contrast, only a single wAna variant was found among the seven integrated wAna genomes examined in lines from Africa, south Asia, and south Pacific islands suggesting that the integration occurred once from a single infectious wAna genome and then spread geographically. Further analysis revealed that for all D. ananassae we examined with the integrated wAna genomes, the majority of the integrated wAna genomic regions is represented in at least two copies suggesting a double integration or single integration followed by an integrated genome duplication. The possible evolutionary mechanism underlying the widespread geographical presence of the duplicate integration of the wAna genome is an intriguing question remaining to be answered.

The Drosophila bag of marbles Gene Interacts Genetically with Wolbachia and Shows Female-Specific Effects of Divergence.

Many reproductive proteins from diverse taxa evolve rapidly and adaptively. These proteins are typically involved in late stages of reproduction such as sperm development and fertilization, and are more often functional in males than females. Surprisingly, many germline stem cell (GSC) regulatory genes, which are essential for the earliest stages of reproduction, also evolve adaptively in Drosophila. One example is the bag of marbles (bam) gene, which is required for GSC differentiation and germline cyst development in females and for regulating mitotic divisions and entry to spermatocyte differentiation in males. Here we show that the extensive divergence of bam between Drosophila melanogaster and D. simulans affects bam function in females but has no apparent effect in males. We further find that infection with Wolbachia pipientis, an endosymbiotic bacterium that can affect host reproduction through various mechanisms, partially suppresses female sterility caused by bam mutations in D. melanogaster and interacts differentially with bam orthologs from D. melanogaster and D. simulans. We propose that the adaptive evolution of bam has been driven at least in part by the long-term interactions between Drosophila species and Wolbachia. More generally, we suggest that microbial infections of the germline may explain the unexpected pattern of evolution of several GSC regulatory genes.

Correction: A nutrient-driven tRNA modification alters translational fidelity and genome-wide protein coding across an animal genus.

[This corrects the article DOI: 10.1371/journal.pbio.1002015.].

Adaptive evolution of genes involved in the regulation of germline stem cells in Drosophila melanogaster and D. simulans.

Population genetic and comparative analyses in diverse taxa have shown that numerous genes involved in reproduction are adaptively evolving. Two genes involved in germline stem cell regulation, bag of marbles (bam) and benign gonial cell neoplasm (bgcn), have been shown previously to experience recurrent, adaptive evolution in both Drosophila melanogaster and D. simulans. Here we report a population genetic survey on eight additional genes involved in germline stem cell regulation in D. melanogaster and D. simulans that reveals all eight of these genes reject a neutral model of evolution in at least one test and one species after correction for multiple testing using a false-discovery rate of 0.05. These genes play diverse roles in the regulation of germline stem cells, suggesting that positive selection in response to several evolutionary pressures may be acting to drive the adaptive evolution of these genes.